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Abstract. Let P be a set of points in R'', and let a ^ 1 be a real number. We define the 
distance between two points p,q £ P as \pq\" , where \pq\ denotes the standard Euclidean 
distance between p and q. We denote the traveling salesman problem under this distance 
function by Tsp(d, a). We design a 5- approximation algorithm for Tsp(2,2) and generalize 
this result to obtain an approximation factor of 3"~^ -f \/6°/3 for d = 2 and all q ^ 2. 

We also study the variant Rev-TSP of the problem where the traveling salesman is 
allowed to revisit points. We present a polynomial-time approximation scheme for Rev- 
Tsp(2,a) with a ^ 2, and we show that Rev-Tsp(d, a) is APX-hard if d ^ 3 and a > 1. 
The APX-hardness proof carries over to Tsp(d, a) for the same parameter ranges. 



1. Introduction 

Motivated by a power-assignment problem in wireless networks (see below for a short 
discussion of this application) Funke et al. [12] studied the following special case Tsp((i, a) 
of the Traveling Salesman Problem (Tsp) which is specified by an integer d ^ 2 and a real 
number a > 0. The cities are n points in d-dimensional space M"^, and the distance between 
two points p and q is \pq\", where \pq\ denotes the standard Euclidean distance between p 
and q. 

• The objective in problem Tsp(d, a) is to find a shortest tour (under distances | • |") 
that visits every city exactly once. 
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• In the closely related problem Rev-Tsp((i, q), the objective is to find a shortest tour 
that visits every city at least once; thus the salesman is allowed to revisit cities. 
Note that Tsp(2,l) is the classical two-dimensional Euclidean TSP and that Tsp((i, oo) is 
the so-called bottleneck TsP in W^, where the goal is to find a tour whose longest edge has 
minimum length. We are, however, mainly interested in the case where a is some small 
constant, and we will not touch the case a = oo. 

Similarities and differences to the classical Euclidean TSP. The classical Euclidean Tsp is 
NP-hard even in two dimensions, but it is relatively easy to approximate. In particular, it 
admits a polynomial-time approximation scheme: Given a parameter e > and a set of n 
points in d-dimensional Euclidean space, one can find in 2^^) + ((i/e)0('^)nlog n time a 
tour whose length is at most 1 + e times the optimal length [23] . 

A crucial property of the Euclidean Tsp is that the underlying Euclidean distances 
satisfy the triangle inequality. The triangle inequality implies that no reasonable salesman 
would ever revisit the same city: Instead of returning to a city, it is always cheaper to skip 
the city and to travel directly to the successor city. All positive approximation results for 
the Euclidean Tsp rely heavily on the triangle inequality. In strong contrast to this, for 
exponents a > 1 the distance function | • |° does not satisfy the triangle inequality. Thus 
the combinatorial structure of the problem changes significantly — for example, revisits may 
suddenly become helpful — and the existing approximation algorithms for Euclidean Tsp 
cannot be applied. 

Another nice property of the classical Euclidean problem Tsp(2, 1) is that, sloppily 
speaking, instances with many cities have long optimal tours. Consider for instance a set P 
of n points in the unit square. Then there exists a tour whose Euclidean length is bounded 
by 0{^/n) |15j . This bound is essentially tight since there are point sets for which every 
tour has Euclidean length ^}{^/n). Interestingly, these results do not carry over to Tsp(2, 2) 
with squared Euclidean distances. Problem #124 in the book by Bollobas [8j shows that 
there always exists a tour for P such that the sum of the squared Euclidean distances is 
bounded by 4, and that this bound of 4 is best possible. Since, as a rule of thumb, large 
objective values are easier to approximate than small objective values, this already indicates 
a substantial difference in the approximability behaviors of Tsp(2,l) and Tsp(2,2). 

Previous work and our results. Funke et al. note that the distance function | • |" satisfies 
the so-called r-relaxed triangle inequality with parameter r = 2"~^ (see Section [2] for a 
definition). The classical TSP under the r-relaxed triangle inequality has been extensively 
studied [21 [3l [U [7] , and all the corresponding machinery from the literature can be applied 
directly to Tsp(ci, a). For instance, Andreae [6] derives a (r^ + r)-approximation for the 
classical Tsp under the r-relaxed triangle inequality (A,-- Tsp, for short). This result trans- 
lates into a (4"~^ -|- 2"^-'^)-approximation for Tsp(- ,a). For r > 3, it is better to apply 
Bender and Chekuri's 4r-approximation [2j for At— Tsp, which yields a 2""''^-approximation 
for Tsp( • , a). Funke et al. derive a (2 • 3°~^)-approximation algorithm for Tsp( • , a), which 
for the range 2 < a < log3/2 3 ~ 2.71 is better than applying the known results [6l |2]. 
The best result for a < 2 is obtained by Bockenhauer et al. [7] whose Christofides-based 
(3r^/2)-approximation for A,-- TsP yields a (3 • 2^"~^)-approximation for Tsp(- ,a). 

We will demonstrate in Section [2] that essentially every variant of the orig inal T^- 
algorithm by Andreae and Bandelt [3] already gives a (2-3"~^)-approximation for Tsp(d, a). 
The bottom-line of all this, and the actual starting point of our paper, is that the machinery 
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around the r-relaxed triangle inequality only yields a bound of roug hly 2 • 3°-^ This 
raises the following questions: How much can geometry help us in getting even better 
approximation ratios? Can we beat the 6-approximation for Tsp(2,2) of Funke et al.? We 
answer these questions affirmatively: We develop a new variant of the T'^-algorithm which 
we call the geometric -algorithm. An intricate analysis in Section [3] shows that this yields 
a 5-approximation for Tsp(2,2). We then extend our analysis to Tsp(2,a) with a > 2, 
and thus obtain a (3"~^ + /3)-approximation; see Sectional This new bound is always 
better than the bound 2 • 3"^^ of Funke et al. and of our analysis of the T^-algorithm. 

Finally, in Section [5l we turn our attention to the following two questions: (a) How 
does the approximability of Tsp behave when we make a larger than one? (b) Does al- 
lowing revisits change the complexity or the approximability of the problem? As we know, 
classical Euclidean Tsp (that is, Tsp((i, 1)) is NP-hard [19] and has a polynomial-time ap- 
proximation scheme (PTAS) in any fixed number d of dimensions [1]. On the other hand, 
Rev-Tsp(d, a) has — to the best of our knowledge — not been studied before. Concerning 
question (b), complexity behaves as expected: Rev-Tsp(d, a) is NP-hard for any d ^ 2 and 
any a > 0, and our (straightforward) hardness argument also works for Tsp((i, a). In terms 
of approximability, we show that whereas the two-dimensional problem Rev-Tsp(2, a) still 
has a PTAS for all values a ^ 2, the problem becomes APX-hard for all a > 1 in three 
dimensions. We were surprised that the APX-hardness proof, too, carried over to Tsp(3, a) 
for all a > 1. This inapproximability result stands in strong contrast to the behavior of the 
classical Euclidean Tsp (the case a = 1). 

The connection to wireless networks. Consider a wireless network whose nodes are equipped 
with omni-directional antennas. The nodes are modeled as points in the plane, and every 
node can communicate with all other nodes that are within its transmission radius. The 
power (that is, the energy) needed to achieve a transmission radius of r is roughly propor- 
tional to r° for some real parameter a called the distance-power gradient. Depending on 
environmental conditions, a typically is in the range 2 to 6 \13\ Chapter 1]. The goal is to 
assign powers to the nodes such that the resulting network has certain desirable properties, 
while the overall power consumption is minimized. A widely studied variant has the objec- 
tive to make the resulting network strongly connected [U [TT\ [T6] . Other variants (finding 
broadcast trees; having small hop diameter; etc) have been studied as well. Funke et al. |12j 
suggest that it is useful to have a tour through the network, which can be used to pass a 
virtual token around. The resulting power-assignment problem is Tsp(2,a). 

Another setting related to Tsp(2, a) is the following. Instead of omni-directional anten- 
nas, some wireless networks use directional antennas. This achieves the same transmission 
radius under a smaller energy consumption [171 122j . To model directional antennas, Cara- 
giannis et al. [9] assume that a node can communicate with other nodes in a circular sector 
of a given angle (where the sector's radius is still determined by the power of the node's 
signal). For directional antennas one not only has to assign a power level to each node, but 
also has to decide on the direction in which each node transmits. If the opening angle tends 
to zero and the points are in general position, a strongly connected network becomes a tour. 
Hence, our results on Tsp(2,a) may shed some light on the difficulty of power assignment 
for directional antennas with small opening angles. 
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2. Approximating Tsp( ■, a) 

In this section we lay the basis for our main contribution, a 5-approximation for 
Tsp(2,2) in Section [3l We review known algorithms for a related version of TSP, which 
can be applied to our setting. As it turns out, these algorithms already yield the same 
worst-case bounds as the algorithm that Funke et al. |12) gave recently. 

We recall some definitions. Let 5 be a set, let dist(-, •) : S x S ^ Mj^o be a distance 
function on S, and let r ^ 1. We say that dist(-, •) fulfills the r-relaxed triangle inequality 
if any three elements p,q,r G S satisfy dist(p, r) ^ r • (dist(p, g) + dist(g, r)). Recall that 
we denote by A,— Tsp the Tsp problem on complete graphs whose weight function (when 
viewed as a distance function on the vertices) fulfills the r-relaxed triangle inequality. The 
following lemma, which has been observed by Funke et al. [12], allows us to apply algorithms 
for At-Tsp to our problem. The proof relies on Holder's inequality. 

Lemma 2.1 (|12j). Let a > be a fixed constant. The distance function | • |° : M"' x — >■ 

I^^O) (P) q) ^ IpqI" fulfills the T -relaxed triangle inequality for r = 2"^-*^. 

Andreae and Bandelt [3] gave an approximation algorithm for A,-- Tsp. Their T^- 
algorithm is an adaptation of the well-known double-spanning-tree heuristic for Tsp. This 
heuristic finds a minimum spanning tree (MST) in the given graph G, doubles all edges, 
finds an Euler tour in the resulting multigraph, and finally constructs a Hamiltonian cycle 
from the Euler tour by skipping all nodes that have already been visited. The weight of 
the MST is a lower bound for the length of a TsP-tour since removing any edge from a 
tour yields a spanning tree whose weight is at least the weight of the MST. Note that 
this statement holds for arbitrary weight functions. If the triangle inequality holds, the 
heuristic yields a 2-approximation since then skipping over visited nodes never increases 
the length of the tour, which initially equals twice the weight of the MST. For the weight 
function | • |", however, the heuristic can perform arbitrarily badly — consider a sequence of 
n equally-spaced points on a line. 

The T^-algorithm of Andreae and Bandelt also creates a Hamiltonian tour by short- 
cutting the MST, but their algorithm never skips more than two consecutive nodes. It is 
never necessary to skip more than two consecutive nodes because the cube of a tree T is 
always Hamiltonian by a result of Sekanina [2l]. Recall that the cube of a graph G contains 
an edge uv if there is a path from u to v in G that uses at most three edges. The proof of 
Sekanina is constructive; Andreae and Bandelt use it to construct a tour in MST'^. 

The recursive procedure of Sekanina [23] to obtain a Hamiltonian cycle in intuitively 
works as illustrated in Fig. [U for the pseudo-code, see Algorithm [TJ The algorithm is 
applied to a tree T and an edge e = uiU2 of T. Removing the edge e splits the tree into 
two components Ti and T2. In each component Tj (z = 1,2), the algorithm selects an 
arbitrary edge Cj = UiWi incident to ut and recursively computes a Hamiltonian cycle of Tj 
that includes the edge €{. The algorithm returns a Hamiltonian cycle of T that includes e. 
The cycle consists of the cycles in Ti and T2 without the edges ei and 62, respectively. The 
two resulting paths are stitched together with the help of e and the new edge wiW2- 

Note that different choices of the edge Cj in line [5] give rise to different versions of the 
algorithm. The standard T^-algorithm takes an arbitrary such edge, while Andreae's refined 
version [2] makes a specific choice, which gives a better result. (In the next section we will 
choose Cj based on the local geometry of the MST, which will lead to an improved result 
for our problem.) Andreae's tour in MST^ has weight at most (r^ + r) times the weight of 
the MST, which is worst-case optimal Combining his result with Lemma [2.11 yields that 
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Algorithm 1: CYCLElNCuBE(r, e = ^1^2) 




Ml 1*2 .. 



1 for i ^ 1 to 2 do 

2 Ti -^r- component of T — e that contains Ui 

3 if \Ti\ = 1 then Pi ^ 0; -Wj ^ Ui 

4 else 

5 pick an edge = tijUij incident to Ui in Tj 

6 if \Ti\ = 2 then Ui ^ Cj 

7 else Ilj ^ CYCLElNCuBE(Tj, Cj) — 



Figure 1 



Recursively find- 
ing a Hamiltonian 
cycle in the cube 
of the tree T. 



8 return Hi + e + 112 + W1W2 



the refined T'^-algorithm is a (4"^^ + 2"^^)-approximation for Tsp( • , a). We now improve 
on this with the help of a simple argument. We will frequently use the following definition. 
Let T be a tree and let vq, . . . ,Vk be a simple path in T. Then we call VQVk a k-shortcut 
of T. We say that a shortcut vw uses an edge e if e lies on the path connecting v and w 
in T. It is not hard to see that the weight of a A;-shortcut can be bounded as follows. 

Lemma 2.2. Let a ^ 1 and let e be a k-shortcut using edges ei,...,ek- Then |e|" ^ 



Given a tree T, the tour constructed by the T^-algorithm consists of edges of T and 2- 
and 3-shortcuts that use edges of T. Note that in this tour each edge of T is used exactly 
twice. Thus, for a ^ 2, the original T^-algorithm does actually better than the bound we 
obtained above for the refined T'^-algorithm. 

Corollary 2.3. Every version of the -algorithm is a {2-3°'~^)- approximation for Tsp( ■ ,a). 

Note that our improved analysis of the T^-algorithm yields the same result as the 
algorithm of Funke et al. [12J. 

Bender and Chekuri [6] designed a 4r-approximation for A,— Tsp using a different lower 
bound: the optimal Tsp tour is a biconnected subgraph of the original graph. The weight 
of the optimal Tsp tour is at least that of the minimum-weight biconnected subgraph. 
The latter is NP-hard to compute [lOj, but can be approximated within a factor of 2 [21j. 
Moreover, the square of a biconnected subgraph is always Hamiltonian. Thus using only 
edges of the biconnected subgraph and iwo-shortcuts yields a 4r-approximation for A^- 
Tsp. Combining the result of Bender and Chekuri with Lemma l2 . II immediately yields the 
following result, which is better than Corollary [2^ for a > log3/2 3 w 2.71. 

Corollary 2.4. The algorithm of Bender and Chekuri is a l"''^^ -approximation for Tsp( ■ ,a). 

3. A 5-Approximation for TSP(2,2) 

In the previous section we have used graph-theoretic arguments to determine the per- 
formance of the T^-algorithm. By Corollary 12.31 the T^-algorithm yields a 6-approximation 
for a = 2, independently of the dimension of the underlying Euclidean space. We now 
define what we call the geometric T'^-algorithm and show that it yields a 5-approximation 
for Tsp(2, 2). The geometric T^-algorithm simply chooses in line [5] of Algorithm [T] the edge 
Cj that makes the smallest angle with the edge e. 



a 
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(a) a and c lie on the same side of the 
line through b 




(b) a and c lie on different sides of the line 
through b 



Figure 2: Two cases for computing the length of the 3-shortcut s{a,b,c) 



The idea behind taking advantage of geometry is as follows. In Corollary 12.31 we have 
exploited the fact that each edge is used in two 3)-shortcuts. The weight of a 3-shortcut 
is maximum if the corresponding points lie on a line. For the case of the Euclidean MST it 
is well-known that edges make an angle of at least 7r/3 if they share an endpoint. The same 
proof also works for the MST w.r.t. | • |". This guarantees that in line[5]of Algorithm [H we 
can pick an edge that makes a relatively small angle with e — if the degree of Ui is larger 
than 2. Otherwise, it is easy to see that is used by a 2)- and a 3)-shortcut, which 
is favorable to being used by two 3-shortcuts, see Lemma |2.2[ 

Although the intuition behind our geometric T^-algorithm is clear, its analysis turns 
out to be non-trivial. We start with the following lemma that can be proved with some 
elementary trigonometry. Given two line segments s and t incident to the same point, we 
denote the smaller angle between s and t by Zst and define ipst = — ^st. 

Lemma 3.1. Given a tree T, the 3-shortcut s{a,b,c) that uses the edges a, b,cofT in this 
order has weight 

\s{a,b, c)p = |a|^ -I- -I- |c|^ -I- 2|a||6| cos^/^fea + 2|6||c| cos-^^c + 2|a||c| cos(V'fea + <J • ^fec), 
where 5 = +1 if a and c lie on the same side of the line through b, and 5 = —1 if a and c lie 
on opposite sides. Moreover, |s(a,6,c)p ^ 2|ap -|- |6p -|- 2|cp -|- 2|a||6| cos V'ba + 2|&||c| cos 'i/'ftc- 

Lemma 13.11 (illustrated in Fig. [2]) expresses the weight of a 3-shortcut in terms of the 
lengths of the edges and the angles between them. Now we show that if an edge a is used in 
two 3-shortcuts, two of these angles are related. Note that the T^-algorithm generates the 
two 3-shortcuts that use a in two consecutive recursive calls, see Fig. O The T^-algorithm 
is first applied to edge b and then recursively to edge a. In the recursive call, the shortcut 
s{e,a,d) is generated where d is an edge incident to both a and b. Then the algorithm 
returns from the recursion and generates the 3-shortcut s(a, b, c). Thus a is the middle edge 
in one 3-shortcut and the first or last edge in the other 3-shortcut. We rely on the following. 

Lemma 3.2. If the geometric T'^ -algorithm generates the two 3-shortcuts s{a,b,c} and 
s(e, a, d) in two recursive calls and d is incident to both a and b, then ^i,a ^ (vr — '4)ad)/'^- 

Now we are ready to prove the main result of this section. 

Theorem 3.3. The geometric -algorithm yields a 5- approximation for Tsp(2,2). 

Proof. We express the length of each shortcut s of the T^-tour in terms of the lengths of the 
MST edges that s uses. Changing the perspective, for each MST edge a, we use contrib(a) 
to denote the sum of all terms that contain the factor \a\. The edge a is used in at most 
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two shortcuts. Bounding their lengths yields an upper bound on contrib(a). The sum of all 
contributions relates the length of the T^-tour to that of the MST (w.r.t. | • |°), which in 
turn is a lower bound for the length of an optimal Tsp tour. 

Due to Lemma [2?2l contrib(a) ^ 5|ap if a is used in a (^2)-shortcut on one side and a 
(^3)-shortcut on the other side. So we focus on the case that a is used in two 3-shortcuts, 
see Fig. [3l We rewrite the composite terms in the bound for s{a,b,c) in Lemma |3. II using 
Young's inequality with e, which, given x, y G R and e > 0, states that xy ^ x^/ (2e)+y^e/2. 

Let V be the vertex that is incident to edges a and b. If there are multiple 3-shortcuts 
that use edges that are incident to v then the T'^-algorithm generates these in consecutive 
recursive calls. We renumber the edges incident to v such that the algorithm is first applied 
to vvi, then recursively to VV2 etc. Then there is some i ^ 1 such that b = vvi and 
a = vvi+i because the algorithm is first applied to b and then recursively to a. We define 
■04 = '>Pvvi,vvi+i{= '^ba)- We rewrite the term 2|a| |5| cos ■i/'ba in the bound for |s(a,6,c)p in 
Lemma |3. II as follows. 

2|a||6| cos V'fta = 2\vvi\\vvi+i\cosiii ^ f{\vvi+i\, \vvi\,ipi), (3.1) 
where 

ifV'i^f, 
\vvi\^ + cos^ Vi if V'i < f and 



f{\vVi+l\, \vVi\,llJi 



2 

1 or (z > 1 and ^ f))) 



\vv, 



+ costpi if V'i < f and i > 1 and ipi-i < 



2 ■ 



The second case of inequality ()3.ip follows from Young's inequality with e = 1/ cos ipi and 
the third case from Young's inequality with e = 1. Replacing 2|6||c| cos '06c in the bound for 
|s(a, 6, c)p in Lemma l3.1l is analogous. Together, the two replacements yield the bound 

\s{a,b,c)\^ ^ 2\af + \b\^ + 2\cf + /(|a|, |6U;,a) + /(|c|, I^U^J. (3.2) 

We use (j3.2p to bound the weights of all 3-shortcuts. The weight of the final tour is the 
sum of the weights of all shortcuts. In this sum we can take the two occurrences of an 
edge a = wfj+i together and analyze the contribution of a to the tour. Note that the result 
of (13. 2p is still at most 3(|ap -|- |6p -|- |cp). So if an edge a is used in a (^3)-shortcut on one 
side and a (^ 2)-shortcut on the other side, then we still have that contrib(a) ^ 5|ap. It 
remains to consider the case that a is used in two 3-shortcuts. Let 5(0, b, c) and ^(e, a, d) be 
these 3-shortcuts. The algorithm is first applied to edge b and generates shortcut s{a,b,c), 
where a is the first or the third edge of the shortcut. Then the algorithm is recursively 
applied to edge a and generates shortcut s(e, a, d), where o is the middle edge. Fig. [3] shows 
how the vertices are numbered in this case. 




Figure 3: Two 3-shortcuts that use edge a. 



Figure 4: Illustration of case III. 



246 M. DE BERG, F. VAN NIJNATTEN, R. SITTERS, G. J. WOEGINGER, AND A. WOLFF 



Let (Ta be a function that takes a sum of terms and returns the sum of all terms that 
contain \a\. We derive the following expression for contrib(a). 

contrib(a) = da (weight (s (a, b, c))) + cJa (weight (s(e, a, d))) 

^ aa (2|ap + + 2\cf + f{\a\,\b\,i^ba) + f{\c\, |&Uf,c)) 

+ aa (2|e|2 + |a|2 + 2|d|2 + /(|e|, |aUae) + f{\d\, |aUad)) 
^ 4|ap + aa{f{\vvi+i\, \vvi\,tpi)) + aaif{\vvi+2\, \vvi+i\,ipi+i)) (3.3) 
By definition of / we have to consider three cases in (13. 3p for contrib(a). 
Case I: Tpi ^ 7r/2 or V'j+i ^ vr/2. 

We assume w.l.o.g. that ipi ^ 7r/2. Then we know that /(IfWj+il, {vviljipi)) = and in 
the worst case cra{f{\vvi+2\, \vvi+i\,ipi+i)) ^ |ap. Thus we have that contrib(a) ^ 5|ap. 

Case II: ipi < tt/2 and tpi+i < tt/2 and (i = 1 or (i > 1 and V'i-i ^ '^Z'^))- 
By definition of / we have: 

(^aifi\vVi+l\,\vVi\,'llJi)) = fJa (|vUjp + IWi+lpCOS^ Vi) = japCOS^V'i 

0-a(/(bWi+2|, bWi+lNV'i+l)) = CTa (d^^^i+lP + k^'i+2p) COS V'i+i) = japCOsV'j+l 

Lemma 13.21 states that ipi ^ {tt — V'i+i)/2. We also know that ipi ^ tt hy definition. Thus 
we have 

contrib(a) ^ ^4 + cos^ ^ ^"^"^ + cos Tpi+i^ ^ 5|ap. 

Case III: ipi < tt/2 and V'i+i < "^72 and i > 1 and ipi^i < 7r/2. 

It can be shown that this leads to a contradiction, see Fig. d] (on page 1245]) . 

In cases I and II, the contribution of any edge \a\ to the tour is at most 5|ap. The 
theorem follows by summing up the contributions of all edges. ■ 

When using the MST as a lower bound in the analysis, there is not much room for 
improvement. There are instances of Tsp(2,2) where the T^-algorithm yields a tour whose 
weight is 4^ times that of the MST; see also [TF, Theorem 4.19]. 

4. Approximating Tsp(2,Q!) with a^2 

In this section we generalize the main result of the previous section to a ^ 2. Our new 
bound is always better than the bound 2 • 3°^"^ of Funke et al. [12j, see also Corollarv 12.31 
For Q < 3.41 our bound is better than the bound 2°''^^ that follows from the algorithm of 
Bender and Chekuri see Corollarv 12. 4[ 

Theorem 4.1. The geometric -algorithm yields a (3°~^ + ^/q'^/?))- approximation for 
Tsp(2,a) if a ^2. 

Proof. If an edge a is used in a (^2)-shortcut on one side and a (^3)-shortcut on the other 
side then the total contribution of a to the tour is at most (2"~^ + 3'^~^)|a|" by Lemma [2.21 
So we will focus our analysis again on the case that a is used in two 3-shortcuts. For a = 2 
we can express the weight of a 3-shortcut by Lemma 13.11 and rewrite the composite terms 
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as in inequality (I3.ip . For a > 2 we apply Holder's inequality. 

\sia,b,cr = {\s{a,b,c)\Y' 

^ (2|a|2 + + 2|c|2 + /(|a|, I^U^,) + /(|c|, |6U6c))°^' 

= + P,\b\' + PMT^' (4.1) 

^ 3°/2-i (^/3°/2|a|" + /3°/^|6|'^ + /3,"/2|c|") (4.2) 

We introduced the constants of type /3 to shorten the expression. Note that the last in- 
equality holds only if a > 2. 

In order to bound the contribution of an edge a that is used in two 3-shortcuts we 
follow the proof of Theorem 13.31 Since the assumptions of case III in that proof led to a 
contradiction, it suffices to consider cases I and II. 
Case I: ipi ^ 7r/2 or V'j+i ^ 7r/2. 

contrib(a) ^ 3"/^-^ ((2 + cos^i)"^^ + (2 + cos V'i+i)°/^) \a\'^ 
^ 3"/2-i J^2°/2 +3"/2^ = (^3^-1 + V673) |a|° 
Case II: tl^i < 7r/2 and tpi+i < vr/2 and (z = 1 or (i > 1 and V'j-i ^ ^/2))- 

contrib(a) ^ 3"/2-i (^(2 + cos V'i+i)"/^ + (2 + sin^ V'i+i/2)"/^) la]" 

^ V ' 

Now we use the fact that the function h : [0,27r] — )- M, x 1-^ (2 + cosx)'^ + (2 + sin^x/2)'^ 
attains its maximum value at x = 0. Thus ga also attains its maximum in the range [0, 7r/2) 
in a; = 0. This yields 

contrib(a) ^ 3°/2-i . ^^(0) . |a|" ^ (3"-^ + V673) |a|". 

In both cases we showed that contrib(a) ^ (3"^^ + \/6"/3)|q^|"- The theorem follows for 
a > 2 by summing up the contributions of all edges. The case a = 2 corresponds to 
Theorem 13.31 ■ 



5. The Approximability of TSP and Rev-TSP 

In this section we discuss complexity and approximability of Tsp and its variant Rev- 
Tsp, where the salesman is allowed to revisit the cities. Recall that for any fixed dimension 
d^2, Tsp(d, 1) is NP-hard [l9] and admits a PTAS [Ij. 

Theorem 5.1. Tsp(d, a) and Rev-Tsp{d, a) are i<sp-hard for any d ^ 2 and a > 0. 

Proof. Itai et al. [14] showed that, given n points in the unit grid, it is NP-hard to decide 
whether there is a Tsp tour of Euclidean length n. Thus for both of our problems it is 
NP-hard to distinguish between OPT = n and opt ^ n — 1 + \f2 . m 
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Theorem 5.2. Tsp((i, a) and Rev-Tsp{d, a) are APX-hard for any d ^ 3 and any a > 1. 

Proof. We only discuss the case d = 3 and a = 2 — ah other cases can be settled by slightly 
modified arguments — Tsp. and we only consider Rev-TsP; a similar reduction can be used 
for Tsp. We reduce from {1,2}-TSP, the Tsp on the complete graph where the weight 
of every edge is either 1 or 2; this problem is APX-hard [20]. An instance of {1,2}-TSP 
consists of the complete graph Kn = {Vn,En) with vertex set Vn = {vi,...,Vn}, edge 
set En = {ei, . . . ,em} where m = n{n — l)/2, and edge lengths that are specified by a 
weight function w : En — {1,2}. Given Kn and w, we construct a corresponding instance 
Pn,w C M? of Rev-Tsp(3, 2). 

We start our construction by introducing several auxiliary line segments. For each 
vertex Vi ^ Vn we, define its spine to be the vertical line segment going from point (ni, ni, n) 
to point [ni, ni, nm). For each edge = ViVj G En with i < j, we define two corresponding 
line segments that are parallel to the xy-plane and that are called bones. The first bone 
connects point {ni,ni,nk) on the spine of Vi to the point {nj,ni,nk). The other bone 
connects point {nj,nj,nk) on the spine of vj to the point {nj,ni — Sk,nk), where = 1 if 
w{ek) = 1 and 6k = \/2 if w{ek) = 2. Note that these two bones do not quite touch; they 
are separated by a gap of length Sk . 

In order to get the instance Pn,w of Rev-Tsp(3, 2), we subdivide every single (spine 
or bone) line segment introduced above by a dense, evenly distributed set of points — we 
call these points cities from now on — so that every unit-length piece receives cities. 
The distance between adjacent cities is and so the cost for going from one city to 

an adjacent city is All these cities together form instance Pn,wi and this completes 

our construction. Since we have introduced line segments with a total length of at most 
n ■ n{m — 1) + m ■ 2n{n — 1) < 2n^, the overall number of cities is at most 2n^. 

For 1 ^ i ^ n we call the cities on the spine of Vi and on all bones incident to this spine 
the city cluster of Vi. Traversing all cities within such a city cluster is very cheap; even if 
we visit every city twice, this costs at most 2 • 2n^ /n^^ = 4/n for all cities in all city clusters 
together. In a traveling salesman tour, the only expensive steps occur when the salesman 
jumps from one city cluster to another city cluster. By the above definition of 5^, when 
jumping from bone to bone across the gap corresponding to edge the incurred cost is 
exactly w{ek). Note that jumping from city cluster to city cluster in any other way would 
be much more expensive and would thus not reduce the total cost of the tour. 

Finally, let us show that our reduction is approximation preserving. Fix an e with 
< e < 1. Consider an instance Kn and w of {1,2}-TSP, and assume without loss of 
generality that n > A/e. Consider an optimal tour ttq for this instance. If ttq uses £ ^ 
edges of length 2 and n — i edges of length 1, then it has cost n + i. Given a PTAS for 
Rev-Tsp, we show how to compute in polynomial time a tour of cost at most (1 + e)(n + i) 
for Kn and w. 

First note that the tour ttq can be transformed into a tour tti through Pn.w that makes I 
jumps of cost 2 and n — i jumps of cost 1. That tour vri costs at most n + ^ + 4/n. Using 
our hypothetical PTAS for Rev-TsP, we can compute for any e' > in polynomial time 
a tour 1T2 through Pn^w of cost at most (1 + e')copt, where Copt is the cost of an optimal 
Rev-Tsp tour. The existence of vri yields Copt ^ n+£+A/n. The tour 7r2 can be transformed 
into a tour vrs through Kn- Just map the jumps of 7r2 to the corresponding edges of Kn. 
Since this mapping cannot increase the cost, tour vrs costs at most (1 -|- e')(n + £ + 4/n). 
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Choosing e' = e/2 and using 4/n<e<l, we can bound the cost of vrs from above by 

(l + |)(n + ^) + (l + |). = (l + i)(n + £) + |(2 + e) = (l + e)(n + ^) 

as desired. Like 7r2, the tour 713 may visit vertices more than once. This can be fixed by 
greedily introducing shortcuts. The shortcuts do not increase the cost of the tour since the 
weight function w (triviahy) fulfihs the triangle inequality. ■ 

Theorem 5.3. There exists a PTAS for i?ew-Tsp(2, a) for any a ^ 2. 

Proof. Given a set P of points in the plane, consider the Gabriel graph Gp that has a vertex 
for each point in P. There is an edge between points p and g, if the open disk with diameter 
pq is empty, in other words, if for all points r G P \ {p, q}, the angle Zprq is at most tt/2. 
The weight of the edge is \pq\°'- Note that |pr|" + |rg|° ^ \pq\°' if Zprq is at least tt/2. 
Therefore, there is an optimal Tsp tour with revisits through P that only uses the edges 
of Gp: Indeed, if a tour uses an edge pq for which there is a point r with Zprq ^ vr/2, 
then replacing pq by pr and rq would shorten the tour. Such a replacement is feasible since 
revisiting city r is allowed. The Gabriel graph is planar. Hence we end up with an instance 
of the Tsp on weighted planar graphs, for which a PTAS is known [5j. ■ 

Recall that a quasi-PTAS is an approximation scheme with running time nP°^y^°§", 
where n is the size of the input. The following result follows immediately from the facts 
that (a) the metric | • |" has bounded doubling dimension and (b) Tsp on metrics of bounded 
doubling dimension admits a quasi-PTAS [25| . 

Theorem 5.4. There exists a quasi-PTAS for Rev-Tsp{d, a) for any a G (0, 1] and d ^ 1. 
6. Conclusions 

In order to construct considerably better approximation algorithms for Tsp{d,a), we 
expect that substantially different methods of analysis have to be found. A result of Van 
Nijnatten [181 Theorem 4.19] indicates that there is not much room left for improvement 
as long as we compare to the MST. 

The approximability of Rev-Tsp(2, a) for 1 < a < 2 is an interesting open question. 
We believe that a (quasi)-PTAS may be obtained using the framework of the PTAS for 
weighted planar graph Tsp by Arora et al. [5]. A simple reduction shows that deriving a 
PTAS for our problem is at least as hard as deriving a PTAS for weighted planar graph 
Tsp. Assume we have a PTAS for Rev-Tsp(2, a) for some a > 1. Given a weighted planar 
graph and a planar embedding, we replace each edge by a dense set of points such that 
traversing a subedge basically costs zero. By making one subedge of each edge e longer, we 
can make the cost of that subedge (and thus of e) in Rev-Tsp proportional to the weight 
of e. Then, the costs of the optimal solutions of the two problems will be the same up to an 
arbitrarily small constant factor of 1 + e. Such a reduction is polynomially bounded if all 
weights are polynomially bounded, which can be achieved by a standard rounding scheme. 

A PTAS for Rev-Tsp(2, a) for any a > 1 would be an interesting generalization of the 
existing PTAS's for weighted planar graphs. Ideally, one would have a PTAS with running 
time independent of a since it would contain both Euclidean Tsp and weighted planar 
graph Tsp as special cases. 
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